Discriminating the source of high-energy positrons with AMS-02 
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We study the prospects for discriminating between the dark matter (DM) and pulsar origin of the 
PAMELA positron excess with the Alpha Magnetic Spectrometer AMS-02. We simulate the response 
of AMS-02 to positrons (and electrons) originating from DM annihilations, and determine the pulsar 
parameters (spin-down luminosity, distance and characteristic age) that produce a satisfactory fit 
to the mock AMS-02 data. It turns out that it is always possible to mimic a DM signal with 
pulsars. Although the fit in some cases requires values of spin-down luminosity and characteristic 
age different from those of known pulsars in the ATNF and Fermi-LAT catalogues, these catalogues 
are known to be incomplete, and therefore the pulsar interpretation can hardly be ruled out. We 
also show that if the positron excess is due to a single pulsar, it is always possible to find a DM 
candidate that provides a good fit to the mock AMS-02 data. The discrimination between the two 
scenarios will thus require a better knowledge of the underlying sources, or complementary data. 



I. INTRODUCTION 

The nature of high-energy cosmic ray electrons and 
positrons remains an open problem in modern astro- 
physics and it is currently a matter of intense debate. 
Data accumulated over the years led to the description 
of the electronic component in local cosmic rays as a sin- 
gle power-law with spectral index around 3.4 at multi- 
GeV energies [1], where positrons are much less abun- 
dant than electrons. A standard picture has therefore 
emerged to explain the local flux of cosmic ray electrons 
and positrons in the context of galactic cosmic ray prop- 
agation. In this framework, positrons result from the in- 
clastic scattering of cosmic ray protons and nuclei against 
the interstellar gas in our Galaxy (mainly H and He nu- 
clei), and then diffuse and lose energy before arriving 
at the Earth. Electrons instead are believed to be pri- 
maries, presumably accelerated in supernova remnants 
[2] and injected into the interstellar medium (ISM). In 
the standard picture, the propagation of this injection 
spectrum would lead to the bulk of the local cosmic ray 
electrons. 

Several hints, however, have challenged over the years 
the secondary nature of high-energy cosmic positrons, in 
particular the results of CAPRICE [3] and HEAT [4, 5] 
that indicate a rather large positron fraction at multi- 
GeV energies. In the last few years, a host of observations 
of high-energy electrons and positrons has become avail- 
able, including results from ATIC [6], PPB-BETS [7], 
PAMELA [8, 9], Fermi-LAT [10, 11] andH.E.S.S. [12, 13]. 
PAMELA, for instance, has measured a steeply rising 
positron fraction above ~10 GeV and up to ^100 GeV. 
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Fermi-LAT results, on the other hand, show an electron 
plus positron flux of spectral index ~3.0 above ~20 GeV 
and with a possible hardening at about 300 GeV, while 
H.E.S.S. hints at a cutoff of a few TeV. This so-called 
electron/positron excess is at odds with the standard 
picture described above. In fact, cosmic ray spallation 
on the galactic disk fails to produce enough positrons 
and, more importantly, is incompatible with a positron 
fraction that rises with increasing energies [14, 15]. As 
pointed out in [16], one needs one or more nearby and 
recent sources to accommodate the data. This is be- 
cause ~100 GeV positrons (and electrons) lose energy 
efficiently through inverse Compton scattering and syn- 
chrotron emission, presenting a cooling time of about 
~ 2 x 10 6 yr which translates into kpc-scale diffusion 
distances. Several possible sources were put forward 
many years ago and recently revisited in the light of new 
electron-positron data. 

One of the most popular hypotheses is to invoke dark 
matter (DM) annihilations [17] or decays [18] in the 
galactic halo. Such annihilations or decays would pro- 
duce high-energy electrons and positrons either directly 
or through the decay of secondary particles. Although 
exciting, this interpretation requires non-standard dark 
matter properties - such as high annihilation cross- 
sections [19-24] - that are in tension with other data, 
including radio emission and 7-rays [25-27], or cosmic 
microwave background measurements [28-30]. Moreover, 
the DM interpretation requires annihilation mainly to 
leptonic channels in order to be consistent with the data 
on the antiproton flux, that itself is completely consis- 
tent with what is expected from secondary production 
[17, 31]. Another possible origin for the excess is the 
emission of electrons and positrons from mature pulsars 
[16, 32-35]. Indeed, electrons can be accelerated in the 
magnetosphere of pulsars and, due to the existing mag- 
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netic fields, emit curvature radiation which will gener- 
ate e pairs and subsequently an electromagnetic cas- 
cade. The uncertainties inherent to this scenario are sig- 
nificant, but it is very likely that one or more known 
pulsars contribute non-ncgligibly to the flux of cosmic 
ray electrons and positrons [36] - check Ref. [37] for the 
implications of recent 7-ray observations of supernova 
remnants and pulsar wind nebulae on the positron frac- 
tion (and antiproton-to-proton ratio). A third hypothe- 
sis, put forward in Refs. [38, 39], posits that secondary 
particles are accelerated in the sites where hadronic pri- 
mary cosmic rays are injected. If secondary particles 
(including electrons and positrons, but also antiprotons, 
boron and beryllium) are created in the acceleration site 
by spallation of primaries on the surrounding medium, 
then they will be accelerated themselves to high ener- 
gies. A neat feature of this mechanism is that it can 
be easily tested since, besides a rising positron fraction, 
also rising secondary-to-primary ratios are predicted [40] . 
Still other possibilities exist to accommodate the elec- 
tron/positron excess, see e.g. Section I of Ref. [34] for 
an overview. In the following, we shall consider only the 
DM and the pulsar interpretations. 

The main goal of the present work is to assess whether 
upcoming measurements of cosmic ray electrons and 
positrons will be sufficient to discriminate between dif- 
ferent origins of the electron/positron excess. It is often 
claimed (see e.g. [41]) that DM annihilations directly into 
e + e~ — that produce a sharp spectral cutoff at the mass 
of the DM particle — can be distinguished from a sin- 
gle pulsar spectrum with future data. Here, we start by 
assuming DM annihilations into leptons as the source of 
the cosmic ray lepton excess, anticipate upcoming mea- 
surements and evaluate to what degree one can discard 
the single pulsar hypothesis in that case. The inverse 
problem is also explored: assuming that the origin of 
the excess is a pulsar, and thus that a pulsar-like cut- 
off is detected, we quantify how well one can reject the 
DM hypothesis. Several DM masses are considered and, 
besides direct annihilation into electron-positron pairs, 
democratic annihilation into leptons (33% e + e~, 33% 
33% t + t~) is also considered. In order to study 
the role of known pulsars, we make use of the ATNF cat- 
alogue [42] and the 7-ray pulsars discovered by Fermi- 
LAT [43, 44]. Furthermore, the anisotropy potentially 
produced by individual pulsars is also discussed as a dis- 
criminating tool. 

We model the response of the Alpha Magnetic Spec- 
trometer AMS-02 [45] , scheduled to be installed in the In- 
ternational Space Station in 2011, including both system- 
atic and statistical uncertainties. Future balloon-borne 
and ground-based experiments may also prove useful in 
measuring high-energy electrons and positrons, but we 
shall not consider them here since the associated sys- 
tematics are likely larger than for AMS-02 (due to the 
influence of the atmosphere). 

Several experimental results have been used along 
the work, namely measurements of the electron flux 



(CAPRICE [3], HEAT [46], AMS-01 [47]), the positron 
flux (CAPRICE [3], HEAT [46], AMS-01 [47]), the 
electron plus positron flux (HEAT [46], BETS [48], 
PPB-BETS [7], ATIC [6], H.E.S.S. [12, 13], Fermi-LAT 
[10, 11]), and the positron fraction (CAPRICE [3], HEAT 
[4, 5], AMS-01 [49], PAMELA [9]). 

The paper is organised as follows. In Sec. II, we outline 
the injection spectra of high-energy e^- from DM annihi- 
lation and pulsars, and their propagation through the 
galaxy. In Sec. Ill, we discuss our modelling of AMS-02 
experimental capabilities. In Sec. IV we assess the ca- 
pability of AMS-02 to distinguish between the DM and 
the pulsar hypotheses and, finally, in Sec. V we draw our 
conclusions. 



II. INJECTION AND PROPAGATION OF 
HIGH-ENERGY ELECTRONS AND POSITRONS 

High-energy electrons and positrons in the galactic 
medium are mainly affected by two processes in their 
way to Earth: energy losses and diffusion. Energy losses 
in the multi-GeV range are dominated by inverse Comp- 
ton scattering off cosmic microwave background, optical 
and infrared photons, and synchrotron emission. In the 
Thomson limit, these losses amount globally to b(E) ~ 
b E 2 with & ~ 1-4 x 10" 16 GeV-V 1 . Diffusion instead 
is caused by the galactic magnetic irregularities and it 
is usually parametrized with a homogeneous power-law 
diffusion coefficient D{E) = D (E/GeV) s . Under these 
assumptions and neglecting convection and reaccelera- 
tion, the number density of electrons and positrons per 
unit energy n(x, E, t) is driven by the transport equation 

^-=Q(x,E,t) + D(E)V 2 n+-^[b(E)n] , (1) 

Q being the source term. Usually this equation is solved 
in steady state conditions (dn/dt — 0) and inside a cylin- 
drical diffusive halo of half-thickness L. The local in- 
terstellar flux then follows directly from the solution of 
equation (1) through (p(x Q ,E) — ^n(x Q , E), where x 
is the position of the solar system in galactic coordinates. 

In this section we briefly review the strategies adopted 
to compute the propagated injection spectra of elec- 
trons and positrons from each of the sources studied: 
DM annihilations and pulsar emission. On top of the 
flux generated by each source, one has of course to take 
into account the baseline astrophysical flux described in 
the previous section, where positrons are merely a by- 
product of the spallation of hadronic cosmic rays on 
the disk and electrons derive mainly from a galactic 
primary component. We shall refer to these standard 
yields as "background" since we are interested in the elec- 
tron/positron excess. In order to be roughly compatible 
with Fermi-LAT measurements while not explaining the 
rising positron fraction observed by PAMELA [8, 9], we 
use "model 1" of Ref. [50] as our background reference 
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setup. This model features a common high-energy injec- 
tion index of 2.42 for electrons (above 4 GeV) and nuclei 
(above 9 GV), and the propagation parameters are fixed 
to the following values: L = 4 kpc, D n = 3.6 x 10 28 cm 2 /s 
and S = 0.33. As in [50], we use GALPROP [51, 52] (ver- 
sion 50. lp) to compute the local fluxes of background 
electrons and positrons within the framework just de- 
scribed, and normalise to Fermi-LAT e^- measurements 
at 100 GeV. 

Before arriving at the top of the atmosphere, low- 
energy charged particles such as electrons and positrons 
are modulated by the solar wind. Our analysis focuses 
on energies above 10 GeV so that this effect is minimal. 
Nevertheless, all fluxes are modulated according to the 
force field approximation [53] with a potential (f>p = 550 
MV. 

It is worth stressing that, even though different meth- 
ods are used to compute the propagated fluxes of each 
source and background, the propagation parameters are 
kept fixed to the values mentioned above, and there- 
fore our analysis is performed in a consistent propaga- 
tion framework. However, as extensively discussed in 
the literature, the parameters L, Dq and 5 are not well- 
constrained by present cosmic-ray data, and hence the 
propagation configuration we use is not unique. Adopt- 
ing other propagation parameters would affect differently 
the local fluxes of background, dark matter and pulsar 
electrons and positrons, but our results are representa- 
tive of typical propagation models. Moreover, taking 
into account the uncertainties in the propagation param- 
eters would add flexibility to fit the data and therefore 
strengthen our results. In the future it would be in- 
teresting to study how the AMS-02 ability to discrim- 
inate high-energy electron-positron sources depends on 
propagation details, perhaps using already the cosmic- 
ray nuclei measurements that will hopefully be provided 
by AMS-02 itself. Along these lines, another matter that 
deserves further investigation is the impact of the hard- 
ening of the proton and helium fluxes recently reported 
by CREAM [54] and ATIC [55] at - TcV/n energies (see 
also [56-58]). 

A. Dark matter annihilations 

If dark matter is composed of annihilating Weakly In- 
teracting Massive Particles (WIMPs), the galactic halo 
should be regarded as a continuous, nearby source of elec- 
trons and positrons. WIMPs of mass m DM and total av- 
eraged annihilation cross-section times velocity (a ann v) 
give rise to a constant source term 

Q^E) = \( p J^)\. ann v)Y J BR k %{E) , 
* \ rriDM / ^ o-hj 

(2) 

where Pdm is the density of dark matter, k runs over the 
annihilation final states, BRk is the branching ratio of 
final state k and is the electron or positron energy 



spectrum per annihilation into final state k. The factor 
1/2 in equation (2) is valid for Majorana self-annihilating 
fcrmions. 

The distribution of dark matter in the Milky Way is 
poorly constrained by observations, but one can learn a 
great deal with the help of numerical simulations. In 
particular, one strong prediction of pure cold dark mat- 
ter simulations is a large galactic population of virialised 
clumps in addition to a smooth halo component. Here, 
the DM density p DM (x) is modelled according to the 
high-resolution dark matter only simulation Via Lactea II 
[59] , where the smooth+clumpy distribution is well fitted 
by a Navarro-Frenk- White (NFW) profile and the abun- 
dance of subhalos of mass M is proportional to M -2 . 
Within this setup, the spherically averaged local smooth 
DM density is p & = 0.41 GcV/cm 3 [60, 61], and we ex- 
trapolate the results of the simulation down to a minimal 
subhalo mass of 10~ 6 M Q which is a fiducial value for 
WIMPs [62] even though it may vary by several orders 
of magnitude [63] . Further technical details can be found 
elsewhere [60, 61]. 

Now, the transport equation (1) with source term (2) 
can be solved semi-analytically. We use the formulas de- 
rived in Refs. [64, 65] to compute the local flux of elec- 
trons and positrons created by DM annihilations in the 
smooth and clumpy components. Notice that the so- 
called substructure boost factor — which quantifies the 
enhancement in the annihilation flux due to the pres- 
ence of DM clumps with respect to the smooth only 
case — is not a constant, but rather an energy- and 
particle-dependent function [64]. Also, unlike sometimes 
assumed, it amounts to small rather than large enhance- 
ments. 

We wish to point out that the distribution of DM parti- 
cles in our Galaxy is uncertain, for what concerns the ra- 
dial profile, the substructure population as well as the lo- 
cal DM density. However, such uncertainties are not cru- 
cial for our work and do not change our conclusions since 
different DM distributions or local densities would simply 
correspond to different normalisations (a ann v) and thus 
would produce essentially the same spectral features (un- 
less the unlikely case of a nearby, massive clump is con- 
sidered) . 

Two annihilation channels will be studied in the 
present paper: direct annihilation to e+e~, and demo- 
cratic annihilation to charged leptons (i.e. BR e + e - = 
BRii+ii- = BR T + T - = 1/3). While in the former case 
the injection spectrum of electrons and positrons is sim- 
ply a Dirac delta at E — m DM , the latter model injects 
pairs with a broad range of energies up to ttldm- In- 
deed, muons decay almost entirely to electrons and taus 
decay either to electrons, or to muons and hadronic parti- 
cles that subsequently generate more electrons. In order 
to model these chains and find the e ± injection spectra 
due to annihilations into p + p~ and t + t~, for each value 
of the DM mass we have generated one million annihila- 
tion events using PYTHIA [66-68] and used the results 
to reconstruct the electron/positron spectrum for a sin- 
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gle annihilation. We have repeated this procedure for 
26 values of the DM mass with constant linear spacing 
Am DM = 20 GeV between 100 and 600 GeV. Single an- 
nihilation spectra for intermediate values of the DM mass 
are obtained by interpolation. 



Other annihilation channels such as light quarks or 
gauge bosons are not considered because they lead to less 
characteristic features in the electron-positron spectrum 
and are thus more challenging to discriminate. More- 
over, when normalized to the PAMELA positron frac- 
tion, models with annihilation to light quarks or gauge 
bosons usually produce an unacceptably large antiproton 
flux, at variance with observations [17, 31]. Note that the 
annihilation modes under study, i.e. direct annihilation 
to leptons only, are difficult to realise in the context of 
minimal supersymmetrical models, but are typical of the 
so-called leptophilic DM models [69] , where tree- level DM 
annihilations to states other than leptons are forbidden 
by an ad hoc symmetry (see also [70] for a supersymmet- 
ric implementation of leptophilic DM) . Models of Kaluza- 
Klein (KK) DM also preferentially annihilate to charged 
leptons, in a democratic way (typically 20% of time for 
each family), and despite the significant branching ratio 
to hadronic channels (~ 35%), can satisfy the constraints 
on the antiproton flux [71]. Thus, even if we do not refer 
to any particular particle physics model in the following, 
but rather choose a more phenomenological approach, 
our results concerning the "democratic direct annihila- 
tion" case can be considered as somewhat representative 
of what one would obtain in the context of leptophilic 
and KK DM models, once the differences in the relative 
branching ratios to leptons are taken into account. On 
the other hand, another class of models, inspired by the 
exciting DM scenario [72, 73], are those in which a new, 
light force carrier is introduced [21]. In these models, the 
DM particle does not couple directly to Standard Model 
states, but instead annihilates to the new state, which in 
turn decays to Standard Model particles. If the new state 
is light enough, the decay to hadrons is kinematically for- 
bidden. Thus the lack of hadrons is determined by the 
spectrum of the states and not by some kind of symmetry. 
However, the fact that the production of leptons happens 
through the annihilation to an intermediate state, makes 
the shape of the final e + /e~ spectrum qualitatively dif- 
ferent from the case in which the DM annihilates directly 
to leptons (as in leptophilic and KK models). For this 
reason, care should be taken in interpreting our results 
in relation to this class of models. 



In this framework, for each annihilation model there 
are thus two free parameters in our analysis: the DM 
mass itidm and the annihilation cross-section (a ann v) 
(the latter entering basically as a normalisation factor). 
We shall use both in a model-independent, phenomeno- 
logical manner in the remainder of the work. 



B. Pulsars 

Pulsars are highly magnetised, rotating neutron stars. 
Their main observational feature is the emission of pulsed 
and directional electromagnetic radiation (from radio 
photons to 7-rays) which suggests that the magnetic 
and rotation axes are misaligned. Unlike dark matter, 
pulsars are known sources of high-energy electrons and 
positrons, being thus a natural class of candidates to 
explain the cosmic ray lepton excess as pursued for ex- 
ample in Refs. [33, 34, 50]. In fact, the magnetosphere 
of a pulsar can easily host potential gaps in excess of 
10 12 V which accelerate primary electrons to TeV ener- 
gies or above. These electrons quickly emit synchrotron 
and/or curvature radiation in the strong magnetic field, 
or upscatter background (radio, microwave, infrared or 
X-ray) photons to TeV-scale energies. Such 7-rays in 
turn are likely to interact with existing low energy pho- 
tons or the magnetic fields and create high-energy 
pairs that subsequently generate an electromagnetic cas- 
cade. Further acceleration may also occur in the nebula 
or remnant surrounding the pulsar. Consequently, along 
with 7-rays (that have been detected), pulsars are very 
credible sources of high-energy cosmic ray electrons and 
positrons. 

The rotation frequency of a pulsar Q = lix / P decreases 
with time so that its rotational energy is dissipated at 
a rate E = f t (|/0 2 ) = Iflti, where I = fM^ 2 . is 
the moment of inertia of a spherical pulsar with mass 
and radius Magnetic braking, that accounts for 
the energy lost by magnetic dipole emission, is certainly 
a mechanism contributing to this spin-down behaviour 
even though there might be others. This process amounts 
to an energy loss rate [2] E mag = -(8TTfl 4 RlBf)/(3c 3 no), 
where B* is the magnetic field at the star surface. In the 
simple case where magnetic braking dominates the pulsar 
spin-down luminosity, one can write E ~ E mag that leads 
to 

/ t\- 1/2 
Q(t) = ft (l + -J , (3) 

f2o being the initial (t = 0) rotational frequency and To = 
(3/^ c 3 )/(167ri?6B 2 fig). Notice that t is essentially the 
luminosity decay time since E = 7f2f2 cx (1 + t/-r )~ 2 . 
Therefore, very old pulsars (t To) — whose early 
e emission has already diffused and diluted — are not 
very likely to contribute to the bulk of the present local 
flux. Put another way, local high-energy electrons and 
positrons must have been produced recently as argued 
in Section I. On the other hand, the electron-positron 
pairs produced by young pulsars need to diffuse through 
the pulsar nebula or the surrounding supernova remnant 
before reaching the ISM, which takes to <~ 10 4 — 10 5 yr. 
The two opposite effects make mature pulsars — of ages 
around 10 5 yr — the dominant source of a pulsar-induced 
e 1 * 1 galactic flux, and hence we shall focus on this type of 
pulsars in the present work. Let us notice at this point 
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that the escape of e^- pairs into the ISM is not fully un- 
derstood yet - an important step along this direction has 
been taken in [74]. 

The characteristic age of a pulsar is obtained by inte- 
grating E ~ E mag under the assumption that the initial 
rotational frequency Q is very large (Clo ^> 0). This 
results in the well-known expression t c h = — f2/(2f2) = 
P/(2P). The actual age of the pulsar, tpsn, may however 
differ from t c h, the discrepancy being of order 10 4 — 10 5 
yr [34] . We shall disregard such discrepancy and identify 
tpsn with t ch . 

A major ingredient to compute the yield of electron- 
positron pairs injected by pulsars is the electronic energy 

output E e ±(t c h) — rje± Jo* Ch 1^1' ^ e± being the fraction 
of rotational energy transferred to electrons and positrons 
(we assume rj e + = i] e - = r/ e ± /2). Using equation (3) and 
assuming t c h » t (roughly valid for mature pulsars), 
one obtains 



arbitrary energy losses and injection spectrum. Applying 
that result with the source term (5), the local density of 
electrons and positrons is found to be 



n e ±(d,E,t ch ) = 



Q e± (E')b(E') 
b{E)^/2 r l f{E ,t ch ) { 



where E'(E, t ch ) = E/(l - E/E rnax (t ch )) is the ini- 
tial energy of a particle detected at energy E assuming 
b(E) = b E 2 , E max (t ch ) = 1/ (bot c h) is the maximum en- 
ergy allowed by losses, d is the distance to the considered 
pulsar, and 



J E 



ri if (E,t) - 4 / dE D(E)/b(E) 



AD(E)tE max (t) (-,_(-,_ E 
(l-S)E { { E max (t) 



i-s\ 



E e ±(t c h) ^ Ve±^l ^ Ve± 



\E\-^ 

TO 



(4) 



where in the last step we have used the approximate 
behaviour of equation (3) for t c h 3> To, $l(t c h) — 
(tch/ To) ■ We have assumed here that the output 
in electron-positron pairs is proportional to the rotational 
energy loss. This needs not be the case; for alternative 
scenarios see e.g. Ref. [34]. 

As for the e ± injection spectrum we use a rather 
phenomenological approach by considering a power-law 
with an exponential cutoff at high energies, Q e ± (E) = 
Q® ± (E/GeV)~ r ex.p(—E/E cut ). Since electron-positron 
pairs and 7-rays are produced by the same physical pro- 
cess, the spectral index T is probably correlated to the 
spectral indices seen in 7-ray observations of pulsars, or 
other multiwavclength measurements, which suggest 1 < 
r < 2. The cutoff energy E cut instead is usually placed 
above the TeV. Finally, the normalisation Q° e± is fixed 
by the output (4) through E e ±(t ch ) = J~ dE E Q e ±(E). 

We now turn to the propagation of the injection spec- 
trum just described. Pulsars may be modelled as burst- 
ing, point-like sources of e ± pairs [32-34, 36, 50]. In- 
deed, the injection region is much smaller than the typi- 
cal propagation distance covered by high-energy electrons 
and positrons, and the emission period is much shorter 
than the travel time to Earth. Thus, the source term 
reads 

Q e ±(x,E,t) = Q e ±(E)6(x-r )6(t-to) , (5) 

in which r is the position of the pulsar and we shall con- 
sider times t — t c h ^> h as natural for mature pulsars. 
In addition, the problem assumes spherical symmetry if 
one considers local sources, namely at distances smaller 
than the half-thickness of the diffusive halo L. The effect 
of boundaries at z = ±L has been studied in [34] and 
shown to be negligible for the energies and pulsar ages of 
interest if L > 1 — 2 kpc. The spherically symmetric ana- 
lytical solution of equation (1) was derived in Ref. [16] for 



Replacing Q e ± and b in equation (6), 

Q^ ± (g/GcV)- r 
n e ± (d, E, t ch ) = ~- o I 1 - 



E 



x cxp 



n 3/2 rdif(E,t ch ) V E max (t ch ) 
I Ej E C ut d 2 



r-2 



1 - E/E max (t c h) r dif(E,t c h) 



(7) 



for E < E maxi and otherwise. Notice that n e - = 
n e + — n e ± /2 since we are assuming rj e + = t] e - = r] e ± /2. 
It is worth pointing out that the spectrum in equation 
(7) features a cutoff at min(E max , E cut ). The nature and 
morphology of the cutoff is different in the two cases, 
Emax E cu t ■ For instance, for very large cutoff energies 
E cu t, the maximum energy local electrons and positrons 
can reach is fixed by energy losses and is therefore a func- 
tion of the pulsar age. To illustrate this point we pick 
the three fiducial pulsar setups in Table I that yield the 
positron fractions shown in Figure 1. Notice the different 
cutoff steepness in the each case - we shall return to this 
issue later on. 

Despite the fact that diffusion erases the initial di- 
rection of charged particles, the output of a pulsar still 
creates a dipole anisotropy on the cosmic ray electron- 
positron sky. Given a certain direction in the sky, 
let N max denote the number of detected photons in 
the hemisphere centred in that direction, and N min 
the number of photons in the opposite hemisphere. 
The size of the anisotropy is energy-dependent and 
given by 5 e ± = (N max ~ N min )/(N max + N min ) = 
3D(E)\Vn e ±(E)\/(cnl°J:), n l °£ being the sum of the pul- 
sar signal n e ± and the background. Applying equation 
(7) one readily obtains 

A (HFt \- 6D{E)d ^(d,E,t c h) 

SMd ' E > tch) -crl f (E,tch)nl°t(d,E,t ch ) ' (§) 

The expressions for the electron- or positron-only 
anisotropies are analogous. 
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r 


E cut [GeV] 


\E\ [crg/s] 


tch [yr] 


E max [GcV] d [kpc] 


/ 


PSR1 


1.7 


1000 


10 j& 


5 x 10 b 


453 0.75 


0.9 


PSR2 


1.3 


1000 


4 x 10 34 


12 x 10 5 


189 0.4 


0.9 


PSR3 


1.7 


500 


8 x 10 34 


3 x 10 s 


755 0.2 


0.9 



TABLE I: The pulsar fiducial setups. In all cases to = 10 4 yr and ?7 e ± = 0.4 as explained in the text. / represents the 
background rescaling factor (see Section IV). 
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FIG. 1: The positron fraction produced by the three pul- 
sar fiducial models indicated in Table I. The solid (dashed) 
line shows the total (background) contribution. In order of 
increasing cutoff energies, the plotted curves correspond to 
PSR2, PSR1 and PSR3. 



In order to explore the possible features of the electron- 
positron spectrum generated by mature pulsars we scan 
the parameter space defined by the spin-down luminos- 
ity \E\, the distance d and the characteristic age t ch in 
the ranges \E\ = 10 32 - 10 36 crg/s, d = 0.1 - 5 kpc and 
t ch = 10 4 - 10 7 yr. Following [34, 36], we adopt r = 10 4 
yr and the rather large e ± fractional output ij e ± = 0.4. 
Notice however that t and rj e ± are both degenerate with 
\E\ for mature pulsars — check equation (4) — so that 
our results may be simply rescaled to account for differ- 
ent values of rj e ± for instance. For the spectral injection 
index and the cutoff energy we take the reference values 
r = 1.7 and E cut — 1 TeV, but we also address the effect 
of varying these parameters in the ranges 1.3 < T < 1.9 
and 0.5 < E cut < 10 TeV. Once \E\, d, t ch , r] e ±, r , T 
and E cut are specified, the local flux of high-energy elec- 
trons and positrons produced by a single mature pulsar 
is unambiguously fixed through equation (7). 

Throughout the work we shall use the ATNF catalogue 
[42] , that contains the most comprehensive list of pulsars 
observed in different wavelengths, as well as Fermi-LAT 
7-ray pulsars including both the objects listed in the cat- 
alogue [43] (for the pulsars in the catalogue without dis- 
tance estimate we use the 7-ray distance as determined 



by [36] ) and eight recently discovered ones [44] . This will 
give us insight on the regions in parameter space occu- 
pied by known pulsars and their role in explaining the 
electron/positron excess. Every pulsar catalogue suffers 
from more or less important uncertainties, in particular 
related to the distance and characteristic age estimates. 
Hence, we will show the distance error in Fermi-LAT pul- 
sars and the age uncertainty due to the timing measure- 
ments reported in ATNF catalogue. 

III. MODELLING AMS-02 EXPERIMENTAL 
CAPABILITIES 

In order to study the prospects for discriminating the 
source of the electronic component in cosmic rays, we fo- 
cus on the expected performances of the Alpha Magnetic 
Spectrometer (AMS-02) [45]. Being a large-acceptance 
spectrometer in space, AMS-02 will likely be invalu- 
able in the measurement of cosmic ray nuclei spectra 
and ratios. Moreover, good rejection capabilities will 
enable a precise determination of cosmic ray electrons 
and positrons in the GeV— TeV range. We start by 
studying the energy range 1—300 GeV for both elec- 
trons and positrons [75] , where the energy resolution has 
been determined using a test beam at CERN: ^ = 

^(0.106/ V^/GcV) 2 + (0.0125) 2 [76]. This corresponds 
to - 10.7% (1.4%) energy resolution at E = 1 (300) GeV, 
and conservatively we assume 15 energy bins per decade. 
The recent replacement of the superconducting magnet 
by the permanent one will allow AMS-02 to last more 
than the initial mission duration of 3 years, even though 
with a downgraded performance. In this work we take a 
data-taking period of 18 years. We have, however, veri- 
fied that, for an 1 year period, the projected AMS-02 data 
shown in the following is essentially unchanged since the 
uncertainties are dominated by systematics. 

For a given energy bin of central energy E b and width 
SE, we estimate the detected number of X particles as 

(E'-E b ) 2 

N x = At 5E A x [ dE' 4> x (E r ) 6 , (9) 

where a = AE(Eh)/2, At is the operating time (taken 
to be 18 year as stated in the last paragraph), Ax is the 
geometrical acceptance of the instrument for X particles 
and 4>x is the differential flux. Following [77], we take 
a mean acceptance for electrons and positrons of A e - = 
A e + = 0.045 m 2 sr, valid in the energy range 1—300 GeV. 
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Notice that the Gaussian smearing in the above formula 
is of particular importance for spectral features such as 
cutoffs, while being less relevant for smooth parts of the 
energy spectrum. 

The relative statistical uncertainty inherent to the 
measurement of the flux 4>x in a certain energy bin is 
simply 1 / y/Nx- Systematic errors instead result mainly 
from the misidcntification of other particles as electrons 
and/or positrons. Protons in particular represent the 
major background for the measurement of the electronic 
component. To estimate their number in each energy bin 
we adopt the proton flux measured by BESS [78] (extrap- 
olated where needed) and use the expression in the last 
paragraph without Gaussian smearing and with A p = 0.3 
m 2 sr. For electrons, we adopt a rejection power against 
protons e~ :p = 3x 10 5 [79, 80]. For positrons, the re- 
jection powers against protons and electrons are respec- 
tively e+ : p = 3 x 10 5 and e+ : e" = 10 4 [77]. In a given 
energy bin, the relative systematic uncertainties are thus 



for the electron flux, and 



N p /N e+ N^-/N c 



for 



the positron flux. For the rest of this work we shall add 
in quadrature systematic and statistical uncertainties. 

It is not entirely clear to what extent will AMS-02 mea- 
sure electrons and positrons above a few hundred GeV. 
In addition to the range 1—300 GeV, we also consider 
the window 300—800 GeV with the characteristics de- 
tailed above. Note that this is an optimistic approach 
since the AMS-02 performance will be likely worse above 
a few hundred GeV. 

Finally, we analyse the prospects for detecting a 
dipole anisotropy in the flux of cosmic ray electrons and 
positrons. A nearby source may in fact produce a non- 
negligible anisotropy, in particular at high energies. Ex- 
perimentally, the anisotropy measurement is limited by 
the presence of the (nearly) isotropic electron-positron 
background. Suppose N ± AN is the number of back- 
ground particles (electrons, positrons or both) in a hemi- 
sphere along a certain direction. Here, AJV represents 
a global uncertainty including systematic and statisti- 
cal errors. Then, the minimum detectable anisotropy at 
n sigma is simply do, n <r — nAN/N. Later on we shall 
use the 2a positron and electron plus positron anisotropy 
reaches of AMS-02 in 18 years, modelled according to the 
above formula and the details given in the previous para- 
graphs. We will also apply the 2a electron plus positron 
anisotropy reach of Fermi-LAT (A e ± ~ 1 m 2 sr) after 5 
years and considering statistical errors only, even though 
this is a rather optimistic limit. A dedicated search for 
e ± anisotropies was already performed by the Fermi-LAT 
collaboration with 1 year of data [81] producing interest- 
ing upper limits on 5 e ± . 



IV. DISTINGUISHING DARK MATTER AND 
SINGLE PULSAR SPECTRA 

In this section we quantitatively study the impact of 
future cosmic-ray electron data on distinguishing dark 
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FIG. 2: The electron plus positron spectrum produced by 
DM annihilations into e + e~ (thick solid lines) and democratic 
leptons (thin solid lines). DM masses are fixed to 100 and 300 
GeV, and annihilation cross-sections are chosen arbitrarily for 
plotting purposes. The dashed line indicates the background 
used in this work, corresponding to "model 1" of Ref. [50]. 
The data sets include HEAT [46], BETS [48], PPB-BETS [7], 
ATIC [6], H.E.S.S. 2008 [12], H.E.S.S. 2009 [13] and Fermi- 
LAT [11]. 





m,DM [GeV] 


{(JannV) [cnV/s] 


/ 


DM1 


100 


5.0 x 10~ 2b 


0.97 


DM2 


300 


3.5 x 10" 25 


0.87 


DM3 


500 


9.0 x 10~ 25 


0.83 



TABLE II: The dark matter fiducial setups. Besides mass 
rriDM and annihilation cross-section (a a nnv), also shown is /, 
the background rescaling factor. In all cases direct annihila- 
tion into e + e~ is assumed. 



matter and single pulsars as the source of the lepton ex- 
cess. We start by assuming a dark matter origin and 
assess to what extent one can discard the pulsar hy- 
pothesis. If dark matter particles directly annihilate into 
electron-positron pairs only or democratically into the 
three charged leptons, then a rather abrupt cutoff at the 
DM mass is expected. This is illustrated in Figure 2, 
where the expected e ± background (detailed in Section 
II) is also shown. In the following we shall focus on direct 
annihilations into e since its extreme spectral feature is 
in principle more difficult to mimic with pulsars. We 
consider three phenomenological sets of DM properties 
— summarised in Table II — that feature a rise in the 
positron fraction as seen by PAMELA and a cutoff at 
100, 300 and 500 GeV. The normalisation of the base- 
line e^- flux was rescaled by a factor / (with respect to 
Fermi-LAT data point at 100 GeV) in order to give some 
room for the extra component. All three models produce 
electron plus positron fluxes compatible with Fermi-LAT 
and H.E.S.S. at the 3a level. We use these three models 
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as our fiducial models. 

For each DM fiducial model, "mock" AMS-02 measure- 
ments of the cosmic ray electronic component are gen- 
erated according to the capabilities outlined in section 
III. The resulting positron fraction and positron flux are 
shown in Figure 3. Using these data (we consider only 
energies above 10 GeV to minimise solar modulation de- 
pendence) , the pulsar hypothesis is tested by performing 
a fixed grid scan within the ranges \E\ = 10 32 — 10 36 erg/s, 
d = 0.1 — 5 kpc and t c h = 10 4 — 10 7 yr. In this work we 
are interested in the goodness of fit and not in extracting 
best-fit parameters and the corresponding uncertainties. 
Therefore, we compute the \ 2 f° r each pulsar parameter 
set and draw the contours outside which parameters are 
excluded at 1 — a = 99% or 50% confidence level (using 
as degrees of freedom the number of available energy bins 
above 10 GeV). The two dimensional contours presented 
in the following are obtained by selecting the minimum 
values of x 2 along the hidden dimension. Figure 4 shows 
the results for DM1, T = 1.7 and E cut = 1 TeV, as well 
as ATNF and Fcrmi-LAT pulsars. Also indicated are the 
regions where a single pulsar — plus background — con- 
tributes > 50% to the positron fraction (thick solid red) 
and positron flux (thick dashed red) at ^86 GeV, along 
with the pulsars producing anisotropies visible at 2er to 
Fermi-LAT after 5 years (dotted blue) . The spectra pro- 
duced by the best fit pulsar parameter set is shown by the 
dotted line in the left column of Figure 3 - the distinc- 
tion between DM and pulsars seems virtually impossible 
with AMS-02 mock data. Several comments are in or- 
der here. First, the benchmark DM1 features a cutoff at 
100 GeV which is well inside AMS-02 range and hence 
in this case a precise measurement of the spectral fea- 
ture is anticipated. This is an optimistic scenario where 
the compatible regions in the pulsar parameter space are 
tight, as shown in Figure 4. In particular, the contours 
restrict very effectively the values for the pulsar charac- 
teristic age t c h needed to mimic the DM signal. This is 
because t c h fixes the maximal energy E max ~ l/(bot c h) 
that induces a rather sharp cutoff (in this case where 
E C ut > E max ). The normalisation is instead given by \E\ 
and d: more distant pulsars require larger energy inputs 
to produce the same propagated spectrum. Such a be- 
haviour — illustrated in the top right plot of Figure 4 — 
breaks when d <C r^/ (check equation (7)) in which case 
\E\ alone fixes the normalisation. 

In Figure 5 we show how the confidence regions in the 
plane t c u vs. \E\ are affected by the choice of the spec- 
tral index T, the cutoff energy E cut and the DM fiducial 
model. As it is evident from the left plot of this figure, 
data featuring cutoffs at higher energy select younger pul- 
sars (since E max cx t~ h ), and are compatible with larger 
portions of the pulsar parameter space because of the 
larger experimental uncertainties at high energies. On 
the other hand, changing the spectral index or the cutoff 
energy changes less importantly the allowed regions in 
the pulsar parameter space. 

We summarise in Table III the number of catalogue 



pulsars producing good fits to the mock data and the 
ones expected to contribute non-negligibly to the elec- 
tron/positron anticipated spectrum. Notice that the 
ATNF and Fermi-LAT catalogues have common objects 
(though with different derived properties) and so the fig- 
ures in Table III report the number of different pulsars 
found in each case. The bottom line of this calculation 
is that, from the phenomenological viewpoint, it is pos- 
sible to mimic a DM-like spectrum with single pulsars, 
even in the extreme case of direct DM annihilations into 
electron-positron pairs where a sharp cutoff is present. 
By construction, the mock data shown in Figure 3 are 
better fit by DM than pulsars, but, statistically speak- 
ing, one cannot exclude at a high confidence level the 
regions inside the contours drawn in Figures 4 and 5. 
However, within these regions we find very few or none 
known catalogue pulsars as indicated in Table III. In 
fact, only for the case DM1, T = 1.7, E cut = 500 GeV 
does one find a pulsar setup that can reproduce the 
mock data in a satisfactory manner. This may be in- 
terpreted as an argument against the pulsar hypothe- 
sis, even though catalogues are likely incomplete since 
pulsars emit electromagnetic radiation in a directional 
fashion and thus only a fraction of these objects may 
be observed from Earth - therefore the pulsar hypoth- 
esis cannot be ruled out on this basis. In order to es- 
timate the incompleteness of catalogues, one may com- 
pute the so-called beaming fraction defined as the frac- 
tion of sky each pulsar beam sweeps, F(a ^ (0, tt), A9) = 
± J An dVL=\ [cos (a - AO/2) - cos (a + A0/2)], where 
a is the angle between rotation and magnetic axes, A9/2 
is the half-width of the beam and we have considered 
top-hat beams (see e.g. [82]). Assuming an isotropic 
distribution for a, the mean beaming fraction reads 

(F)(A9) = j*' 2 daF(a,A9)sina = fsm(^). The 
width of the beam A9 is, of course, dependent on the 
emission mechanism. Radio beams, for instance have 
typical half-openings A9/2 = 5.8° (P/s)- 1 / 2 [83] which 
for a fiducial period P <~ 0.1 s gives A9/2 ~ 18° and 
(F r ) <~ 0.24. However, values F r <~ 1 have been found 
by studying different populations of detected 7-ray and 
radio pulsars [82] indicating wide radio emission. As for 
7-rays, beaming fractions in the literature vary in the 
range F g ~ 0.1—1 [84]. Now, if the distribution of 
nearby pulsars is assumed isotropic, the ratio of total to 
observed numbers of objects should be given by 1/F, or 
<~l-4 (1-10) using the above-stated radio (7-ray) beam- 
ing fractions. These figures do not include distance se- 
lection effects, but recall that in the present work we are 
interested in nearby, high-luminosity objects. 

Another difficulty with the single pulsar explanation 
evident from Table III arises from the fact that at least 
a few known pulsars should contribute (together with 
the background) more that 50% of the positron fraction 
and positron flux at a given energy (~ 86 GeV). There- 
fore, in order to accommodate the cutoffs of Figure 3 
with a single pulsar, one also needs to explain why a few 
known pulsars produce less electron-positron pairs than 
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FIG. 3: Mock AMS-02 data (filled blue circles) and corresponding error bars for the DM fiducial setups in Table II. The solid 
(dashed) line shows the total (background) contribution. The top row of frames features the positron flux while the bottom 
one the positron fraction. From left to right, the columns correspond to DM1, DM2 and DM3. In the left column the dotted 
line shows a pulsar spectrum that fits well the AMS-02 mock data. For the positron flux, the data sets include CAPRICE [3 
HEAT [46], AMS-01 [47]. For the positron fraction, the data sets include CAPRICE [3], HEAT 1994-1995 [4], HEAT 2000 5' 
AMS-01 [49] and PAMELA [9]. 
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TABLE III: The number of known pulsars (from ANTF and Fermi- LAT catalogues) producing good fits to the mock AMS-02 
data in Figure 3 for different combinations of P and E cut . jvf 0,99 represent the number of catalogue pulsars found within the 
50%, 99% CL contours, while N2 (N3) is the number of pulsars that contribute individually — plus background — more than 
50% of the positron fraction (flux) at ~86 GeV. 



expected. A caveat to this argument is the fact that a 
rather high efficiency rj e ± = 0.4 is being used. For the 
reference case DM1, V — 1.7, E cut — 1 TeV, we redid the 
calculations with rj e ± = 0.04 and found no known pul- 
sar contributing significantly to the e + flux or fraction, 
which indeed weakens the above-mentioned reasoning. 



It has been proposed in the literature that the pulsar 
origin of the cosmic ray lepton excess could be tested 
by means of anisotropy measurements. Unfortunately, 
however, we find no catalogue pulsars producing e or 
e + anisotropy levels visible to Fermi-LAT in 5 years or 
AMS-02 in 18 years. This conclusion is clear in Fig- 
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FIG. 4: Regions in the pulsar parameter space compatible with the mock AMS-02 data shown in the left column of Figure 3 
corresponding to the fiducial model DM1. The thin inner (outer) solid contour corresponds to an 1 — a — 50% (99%) confidence 
level goodness of fit. The thick solid (dashed) red lines individuate the regions where a single pulsar — plus background — 
contributes >50% of the positron fraction (flux) at ~ 86 GeV. The dotted blue line delimits the 2a e anisotropy reach of 
Fermi-LAT after 5 years. Also shown are the ATNF catalogue [42] (filled dots) and Fermi-LAT pulsars [43, 44] (open circles). 
Notice that not all pulsars visually within the thin contours are actually inside the corresponding three-dimensional confidence 
regions since there is a third, hidden dimension. 



ure 6 where the experimental reaches have been plotted 
along with the maximal anisotropy — attained at energy 
Emax or Emax ~ featured by known pulsars with ages 
t ch = 10 4 — 10 7 yr. Obviously it is possible that a low- 
energy anisotropy is visible while the maximal value goes 
undetected, but we do not consider such case since it 
would not be a strong evidence for singling out a partic- 
ular pulsar. In this figure the references values T = 1.7 
and E cut — 1000 GeV were assumed, but we have also 
checked that the expected anisotropics from known pul- 
sars are below the experimental reaches for all the combi- 
nations of parameters (r, E cut ) listed in Table III. Notice 



that the recent 3a e dipole anisotropy upper limits from 
Fermi-LAT [81] range from - 10~ 2 at 100 GeV to - lO" 1 
at 500 GeV — somewhat in between the dashed and solid 
curves in the left frame of Figure 6. Note as well that 
the anisotropies we find are significantly smaller than in 
other works (see e.g. [85]) because we are using an e^- 
"background" tuned to roughly match Fermi-LAT data 
and we are not letting the catalogue pulsars output \E\ 
vary to explain the lepton excess. 

We now turn to the inverse problem: we assume a 
single pulsar as the source of high-energy electrons and 
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FIG. 5: The 99% CL region compatible with the mock AMS-02 data in Figure 3 for different DM models (left), F (central) 
and E cut (right). 
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FIG. 6: The maximum anisotropy produced by ATNF (filled dots) and Fermi-LAT (open circles) pulsars with t c h = 10 4 — 10 7 
yr as a function of the energy at which such anisotropy is attained. Here, the values F = 1.7, E cut = 1000 GeV and / = 0.9 
(besides to = 10 4 yr and r] e ± = 0.4) were assumed for all catalogue pulsars. The fiducial setups in Table I are shown by red 
circled crosses. In the left (right) frame we show the e (e + ) maximum anisotropy. The AMS-02 18 years (Fermi-LAT 5 years) 

2(7 anisotropy reach is plotted in solid (dashed). Notice that BJj,„ and E^ ax are different from the maximal energy allowed 
by losses, E max . 



positrons, and evaluate the prospects for distinguishing 
this scenario from that of DM annihilations. The three 
fiducial sets of pulsar properties shown in Table I and 
Figure 1 were adopted, each featuring a different cutoff 
sharpness. Applying a procedure very similar to the pre- 
vious case, we generate mock data for the positron frac- 
tion and positron flux in the pulsar scenario and test the 
DM hypothesis. Direct annihilations into e and demo- 
cratic annihilations into leptons are both considered, and 
TTiDM an d (cr ann v) are treated as free parameters. The 
resulting 1 — a = 99%, 50% contours are presented in 
Figure 7. The shaded regions are excluded at 3cr by 
present data on the positron fraction (PAMELA) and/or 
the electron plus positron flux (Fermi-LAT, H.E.S.S.). 



As evident from Figure 7, DM models with direct anni- 
hilations into electron-positron pairs can mimic the spec- 
trum produced by the benchmark PSR2, but not the first 
and third cases that present milder cutoffs. Analogously, 
democratic DM annihilations into charged leptons can- 
not be ruled out as explanation of the spectrum induced 
by PSR1, but do not explain a sharp cutoff as the one 
produced by PSR2. For the sake of clarity we have not 
considered further annihilation channels that may ease 
the fit to mild pulsar-like cutoffs. In any case, for the 
representative pulsar cases, one can always have well-fit 
DM spectra either with annihilations into e ± or demo- 
cratic leptons. Of course, as extensively discussed in the 
literature and already pointed out in Section I, the large 
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log 10 [m D M/GeV] log l0 [m DM /GeV] 

FIG. 7: The regions in the dark matter parameter space compatible with the mock AMS-02 data assuming the three pulsar 
benchmarks of Table I. The inner (outer) contour corresponds to an 1 — a — 50% (99%) confidence level goodness of fit. In the 
left frame, direct annihilations into e + e~ are assumed and the solid contours correspond to PSR2; in this case no compatibility 
is found for PSR1 nor PSR3. In the right frame, democratic annihilations into leptons are assumed and the solid contour 
(1 — a = 99%) corresponds to PSR1; in this case no compatibility is found for PSR2 nor PSR3. In both frames, the shaded 
region is excluded at 3a by PAMELA positron fraction or Fermi-LAT/H.E.S.S. electron plus positron flux, being that the 
portion above the upper solid line is excluded by PAMELA. 



cross-sections highlighted in Figure 7 are excluded by or 
in tension with other data. 

Throughout our work we have assumed (almost) "per- 
fect data", i.e. mock data presenting no fluctuations 
with respect to the "true" observables (but including the 
smearing due to the energy resolution of the instrument). 
This is of course a quite optimistic assumption regarding 
the ability to discriminate different sources and there- 
fore strengthens our conclusions. Nevertheless, for the 
benchmark DM1 we have generated mock data by draw- 
ing the flux in each energy bin from a Gaussian distri- 
bution of mean given by the flux corresponding to Nx 
in equation (9) and standard deviation as explained in 
Section III. In that case, for a single realisation, the true 
model yields a good (but not perfect) fit to the mock data 
with XDM/ n dof — 40.7/58 (p = 95.6%), as opposed to 
Xdm ~ using "perfect data" . On the other hand, fixing 
T = 1.7 and E cut — 1 TeV, we find that the best pulsar 
parameter set provides a chi-square Xp U i/ n dof = 42.1/58 
(p = 94.2%), which cannot be ruled out with any reason- 
able statistical significance. We have also checked that 
the 99% CL contour in the pulsar parameter space ob- 
tained using mock data with fluctuations encompasses a 
region similar (but smaller) than the one shown in Fig- 
ures 4 and 5. In the future, a dedicated study of the per- 
formance of AMS-02 in detecting electrons and positrons 
may allow a more realistic analysis, but we stress here 
that such would reinforce the idea that rejecting the pul- 
sar hypothesis in case of a DM-like spectrum (or vice- 
versa) will be highly challenging. 



V. CONCLUSIONS 

In this work we have studied the capability of fu- 
ture precision measurements to constrain possible ori- 
gins of the cosmic-ray lepton excess, focusing on the 
dark matter and pulsar interpretations. In order to 
model upcoming experimental capabilities, the perfor- 
mance of AMS-02 was estimated in a realistic way, in- 
cluding both systematic and statistical uncertainties. We 
have considered the case in which a cutoff in the elec- 
tron/positron spectrum, produced either by DM annihi- 
lations or by acceleration in a pulsar, is observed, and 
quantified the ability of AMS-02 to reject the wrong hy- 
pothesis. In particular, the first scenario studied was the 
case where the excess is due to DM direct annihilations 
to e pairs, for three different values of the DM mass 
(m DM = 100, 300, 500 GeV). Even if the DM-induced 
spectral cutoff is quite sharp, it turns out that it is al- 
ways possible to find values for the pulsar spin-down lu- 
minosity, distance and characteristic age that produce 
an electron-positron spectrum that would be compati- 
ble, within the experimental uncertainties, with the DM 
one. However, such would require quite specific values 
of the pulsar luminosity and age. In fact, in nearly all 
cases, none of the kwown pulsars within the ATNF cat- 
alogue, nor of the Fermi-LAT pulsars, satisfies these re- 
quirements, although it should be taken into account that 
catalogues are likely incomplete. Then, we have consid- 
ered the symmetric case where the excess is produced by 
a pulsar, and tried to fit the anticipated data with DM 
directly annihilating either to e + /e~ pairs or democrat- 
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ically (i.e. 33% of the time each) to e + /e~, /i + /fi~ and 
t + /t~. Within the framework of each model, DM was 
treated in an as much model-independent way as possi- 
ble, taking the DM mass and annihilation cross section as 
independent parameters; nevertheless, the "democratic" 
case can be seen as somewhat representative of KK and 
leptophilic DM models. We find that the possibility to 
reproduce the pulsar-induced spectrum, as well as the 
volume in the DM parameter space (in other words the 
amount of fine tuning) required to do so critically de- 
pends on the shape of the cutoff, and thus on the charac- 
teristics of the pulsar originating the excess. In any case, 
generically speaking, it is possible to mimic a pulsar-like 
spectrum with DM annihilations, even though the re- 
quired cross-sections are usually in tension with other 
data as extensively discussed in the literature. 

Let us point out that our results were obtained in 
favourable conditions regarding a possible source discrim- 
ination, namely by fixing the propagation setup. Includ- 
ing propagation uncertainties would likely worsen the 
prospects for discrimination. The main conclusion of our 
work is therefore that future e^ data will likely be in- 



sufficient to discriminate between the dark matter and 
the single pulsar interpretations of the cosmic-ray lepton 
excess. One caveat to this statement would be the detec- 
tion of several bumps in the electron-positron spectrum 
at high energies that could be associated to the contribu- 
tion of several nearby pulsars, and that would be difficult 
to mimic with dark matter annihilations or decays (see 
e.g. [86]). In the future, complementary data, or a better 
understanding of both sources and cosmic-ray propaga- 
tion, may contribute to a better discrimination power 
than presented here. 
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